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Abstract 



It is an open question why the temperature of molecular gas in the Galactic center region is higher than 
that of dust. To address this problem, we made simultaneous observations in the NH3 {J,K) = (1,1), (2,2), 
and (3,3) lines of the central molecular zone (CMZ) using the Kagoshima 6 m telescope. The ortho-to-para 
ratio of NH 3 molecules in the CMZ is 1.5-3.5 at most observed area. This ratio is higher than the statistical 
£3 , equilibrium value, and suggests that the formation temperature of NH3 is 11-20 K. This temperature is 

similar to the dust temperature estimated from the submillimeter and infrared continuum. This result 
suggests that the NH3 molecules in the CMZ were produced on dust grains with the currently observed 
temperature (11-20 K), and they were released into the gas phase by supernova shocks or collisions of dust 
particles. The discrepancy between warm molecular gas and cold dust can be explained by the transient 
heating of the interstellar media in the CMZ approximately 10 5 years ago when the NH3 molecules were 
released from the dust. 
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' 1. Introduction {n = 0, 1, 2, •■■) lines. Ortho-NrI 3 and para-NH 3 are kept 
Q\ , for a long time (~ 10 6 yr) because no collisional or fast 
It is one of the biggest questions why the temperature radiative process can transform between these two con- 
of molecular gas in the CMZ is much higher than that of figurations of NH3. The ortho-to-para ratio is expected 
dust. Many molecular line observations have revealed that to be the statistical equilibrium value of 1.0 when NH3 
the temperature of molecular gas is higher than that of molecules are formed in the processes of gas phase reac- 
dust in the Galactic center region (Nagayama et al. 2007, tions under high temperature. The ortho-to-para ratio is 
Paper I and references therein). It is inconsistent with expected to be larger than unity when NH3 molecules ad- 
the standard molecular gas heating model because there sorbed on cold dust grain surfaces are released into the 
interstellar gas is mainly heated up by hot dust. Surveys gas phase with excess desorption energy, which is com- 
in NH3 (J,K) = (1,1) and (2,2) lines have revealed that parable to the energy difference between the ortho and 
this inconsistency is not confined to small regions; in fact, para states (Umemoto et al. 1999). Therefore, the ortho- 
it extends over the entire molecular cloud complex in the to-para ratio indicates the physical conditions in which 
Galactic center, called the CMZ. This implies that a gas NH3 molecules are formed. Moreover, NH3 molecules 
heating process over a hundred parsec scale differs from formed on a dust grain surface colder than 100 K cannot 
that in the galactic disk. be released to interstellar gas without the heating of dust 
To address this problem, the abundance ratio of ortho- (Sandford & Allamandola 1993). NH3 molecules with a 
NH3 to para-NH3, ortho-to-para ratio, is a key. An NH3 low formation temperature indicate that dust should be 
molecule has two different structures. One is called ortho- heated up by some process such as the passage of a shock 
NH3, and the other para-NH3 which arc different in rel- wave (Umemoto et al. 1999). 

ative orientation of the three hydrogen spins with re- In the CMZ, the (3,3) emission is stronger than the 

spect to the nitrogen spin. Ortho-NH3 emits K = 3n (1,1) and (2,2) emissions (Morris et al. 1983). This in- 

(n = 0, 1, 2, •■■) lines while para-NH 3 emits K ^ 3n dicates that ortho-NH 3 is overabundant as compared to 
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para-NH3. However, the ortho-to-para ratio of NH3 over 
the entire CMZ has not yet been obtained. We aim to 
clarify the formation mechanism of molecular gas in the 
CMZ by simultaneous observations in NH3 (1,1), (2,2), 
and (3,3) lines and by estimating the ortho-to-para ratio. 

We present large-scale maps of the CMZ in the NH 3 
(1,1), (2,2), and (3,3) lines. In section 2, we describe our 
observations in detail. The observed data are presented in 
section 3. In section 4, we discuss the physical conditions 
of the CMZ through measurements of the rotational tem- 
perature and the ortho-to-para ratio in NH3 lines. In this 
paper, we assume that the distance to the Galactic center 
is 8.5 kpc and we use the direction on the sky based on 
the galactic coordinates. 

2. Observations 

2.1. Data Obtained Using the Kagoshima 6 m Telescope 

We conducted a large-scale survey using the Kagoshima 
6 m telescope of the National Astronomical Observatory 
of Japan (NAOJ) from September 2005 to March 2007. 
We made simultaneous observations in the NH3 (J, K) = 

(1.1) , (2,2), and (3,3) lines at rest frequencies of 23.694495, 
23.722633, and 23.870129 GHz, respectively. At 23 GHz, 
the telescope beamwidth was 9'.5 and the main beam ef- 
ficiency (?7mb) was 0.54. We used a K-b&nd HEMT am- 
plifier whose system noise temperature was 200-300 K. 
All spectra were obtained using a 8192-channel FX-type 
software-based spectrometer developed by the Kagoshima 
University and the National institute of Information and 
Communications Technology (NiCT). The bandwidth and 
frequency resolution of the spectrometer are 256 MHz 
and 31.25 kHz, respectively. At the NH3 frequencies, 
these correspond to a velocity coverage and velocity res- 
olution of 3200 km s _1 and 0.39 km s^ 1 , respectively. 
We obtained 471 NH 3 (1,1), (2,2) and (3,3) spectra at 
-l.°000 < I < 2.°000 and -0.°375 < b < +0.°375 with a spac- 
ing of 0.°125. The surveyed area corresponds to 450 x 110 
pc. The NH 3 (1,1) and (2,2) data presented in this paper 
were newly observed, although we have already observed 
them, as shown in Paper I. To obtain intensity ratios from 

(2.2) to (1,1) and (3,3) to (1,1) accurately, we simultane- 
ously observed three lines. The data shown in Paper I are 
not included the analysis in this paper. 

All data were obtained by position switching between 
the target positions and reference positions. The reference 
positions are at the Galactic latitude b < — 1°, where no 
NH3 emission was detected. We integrated at least 30 min 
at each point. The relative pointing error was better than 
1', which was verified by the observations of several H2O 
(frequency; 22.235080 GHz) maser sources. 

Data reduction was performed using the UltraSTAR 
package developed by the radio astronomy group at the 
University of Tokyo (Nakajima et al. 2007). To im- 
prove the signal-to-noise ratio, the obtained spectra are 
smoothed to a velocity resolution of 5 or 10 km s _1 . The 
rms noise level after smoothing to 5 km s _1 is typically 
0.030 K in unit of the main beam brightness temperature 
defined by Tmb = T£/ tjmb , where Tj£ is the antenna tem- 



perature calibrated by the chopper wheel method (Kutner 
& Ulich 1981). This sensitivity is better than that in the 
data shown in Paper I by a factor of 2.7. In this paper, the 
intensities are presented in the main beam temperature. 

2.2. Data Obtained Using the Kashima 34 m Telescope 

In order to conduct an investigation of the "0.°9 wing 
feature" , which is a high- velocity wing at I ~ 0.°9 at higher 
resolution (see section 4.5.5), we observed the area at 
0.°700 < I < 0.°940 and -0.°180 < b < +0.°020 with a spac- 
ing of 0.°040. This observation was made using the l'.6 
beam of the Kashima 34 m telescope of NiCT from May 
to November in 2007. The spectra in the NH 3 (1,1), (2,2), 
and (3,3) lines were obtained by the same method as that 
used in the Kagoshima 6 m telescope survey. We inte- 
grated at least 5 min at each point. The rms noise level 
after 5 km s _1 smoothing is typically 0.064 K in Tmb- 

3. Results 

3.1. Profiles 

Direct information about the observed lines can be seen 
in each line profile. Significant NH3 (1,1), (2,2), and (3,3) 
emissions for which the signal-to-noise ratio exceeds 3 were 
detected for 101, 84, and 97 positions out of 157 observed 
positions, respectively. The observed profiles in the (1,1), 
(2,2), and (3,3) lines toward four prominent positions [Sgr 
C (I, b) = (-0P500, -0.°125), Sgr A (0.°125, -0.°125), Sgr 
B (0.°750, -0P125), and l.°3 region (1P250, -0.°125)j are 
shown in Figure 1. The (2,2) line is the weakest, and 
the (3,3) line is the strongest at most positions, especially 
near the velocity of peak intensity. 

The line shapes of the three transitions are similar. 
Although an NH3 line profile has five quadruple hyper- 
fine lines, we cannot separate these components due to 
the large internal motion of the gas in an observed beam. 
The line shapes reflect the violent motion of the gas . 

3.2. Intensity Distributions over the Entire CMZ 

To figure out the NH3 gas distribution on the sky over 
the entire CMZ, the velocity integrated intensity distribu- 
tion is useful for comparison with the maps observed in 
radio and infrared continuum emissions. Figure 2 shows 
the l-b maps of intensities, velocity-integrated in the range 
of wlsr = —200 to 200 km s _1 , covering the NH3 emission. 

The observed area is extended to I > l.°6 from the pre- 
vious survey (Paper I) along the galactic plane. The dis- 
tributions in the (1,1) and (2,2) lines are consistent with 
those shown in Paper I; in addition, fainter peaks and 
diffuse emission were detected because of the improved 
sensitivity. The total integrated intensities in the (1,1), 
(2,2), and (3,3) lines of the entire observed area are 1693, 
1159, and 2120 K km s -1 , respectively. 

The (3,3) line is emitted by the ortho-NH 3 . The (1,1) 
and (2,2) lines are emitted by the para-NH3. Ortho and 
para-NH3 cannot be transformed by collisional or fast ra- 
diative process in the interstellar gas. This implies that 
the ortho- and para-NH3 molecules can be dealt with as 
different molecular species. Therefore, the similarity of 
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the intensity distribution in the (3,3) line with those in 
the (1,1) and (2,2) lines should contain some astrophysi- 
cal information. 

The l-b maps in the NH3 lines appear similar to other 
molecular lines such as those of CO (Sawada et al. 2001) 
and CS (Tsuboi et al. 1999). This similarity is not only in 
the sky distribution but in the l-v domain. Figure 3 shows 
l-v diagrams of three transitions integrated over the ob- 
served latitudes and Figure 4 shows the l-v diagrams at a 
fixed latitude. Even for wlsr = km s _1 , the NH3 lines 
are not affected by foreground absorption; in contrast, the 
effect is significant in the CO line. Due to the difference 
of the critical densities, the CO line traces lower-density 
gas whereas the NH 3 line traces only higher-density gas, 
n(H 2 ) - 10 4 cm" 3 . The NH 3 and CO lines trace the 
different density gas microscopically. Therefore, morpho- 
logical similarity suggests that the flux ratio of core to 
envelope in an observed beam is almost uniform over the 
CMZ. Moreover, the fairly large optical depths in both 
lines (r < 2-4) may reduce the morphological difference, 
because both lines mainly trace the outer envelope of the 
cloud complex. 

3.3. Observational Properties of Individual Clouds 

To investigate the gas properties at various locations in 
the CMZ, wc divide the intensity distribution into several 
small regions and define a cloud as a single peak in the l-v 
and b-v distributions, with the intensities stronger than 
the 3 a level (> 0.1 K in Tmb) in all three lines. Based 
on this definition, we identified 13 clouds. The total inte- 
grated intensity of these 13 clouds is 3447 K km s _1 , which 
corresponds to 69% of the integrated intensity of the en- 
tire observed area. Therefore, the majority of the NH3 gas 
in the CMZ is traced by these clouds. The peak positions, 
apparent sizes, and line widths of these clouds are listed in 
Table 1. We defined their sizes and widths (Al, Ab, and 
Av) as apparent full widths at half-maximum (FWHM) of 
the main-beam temperature without beamsize deconvolu- 
tion. The integrated intensities of the clouds are shown in 
Table 2. 

Four major clouds of the CMZ are seen in the NH3 maps 
(from west to east: the Sgr A 20 km s _1 cloud at I ~ — 0.°1, 
the Sgr A 40 km s" 1 cloud at I ~ 0.°1, the Sgr B cloud at 
I ~ 0.°7, and the l.°3 region cloud at I ~ l.°3). These clouds 
were also detected in the previous observation (Paper I), 
and their locations and extents are found to be the same 
with those of the previous ones. In our new map, we 
clearly detected the NH3 clouds associated with the Sgr 
C HII region at (/, b, v LS r) ^ (-0.°5, -0.°1, -50 km s _1 ) 
and a part of the expanding molecular ring (EMR) at 
(I, b, wlsr) - (-0.°9, -0.°1, 140 km s _1 ). Both had been 
marginally detected in the previous observation. 

3.4- Intensity Ratios 

The intensity ratio of the (2,2) line to the (1,1) line, 
-R(2,2)/(i.i)? is determined by the optical depth, r, and the 
kinetic temperature, 2\, of the gas. With a given optical 
depth, the rotational temperature of the transition, T TOt , 
can be derived from the -R(2,2)/(i,i) an d it gives the gas 



kinetic temperature. Therefore, the changing line ratios 
indicate the differences of the rotational temperature. 

Figure 5(a) shows the histograms of ratio for the in- 
tegrated intensity over the entire observed area with 3 
detection in both the (1,1) and (2,2) lines after smoothing 
to 10 km s _1 . The value of -R(2,2)/(i,i) ranges between 
0.4 and 1.5. The mean value and standard deviation of 
-R(2.2)/(i.i) are derived to be 0.70 and 0.15, respectively. 
The NH3 integrated intensity with -R(2,2)/(i,i) = 0.4-0.8 
and i?( 2 ,2)/(i,i) > °- 8 arc 82% and 18% of the total NH 3 
integrated intensity, respectively. These values are consis- 
tent with the results described in Paper I. 

The intensity ratio of the (3,3) line to the (1,1) line, 
i?(3,3)/(i,i)i is determined by the relative abundance of the 
ortho-to-para NH3, R Q / pi and the kinetic temperature. In 
the case that ortho-NH3 and para-NH3 have the same ki- 
netic temperature, i?(3 j3 )/(i j i) increases with R / p . In the 
case that the ortho-to-para ratio is uniform, R(3 t 3)/n t i) in- 
creases with the kinetic temperature. Wc can estimate the 
ortho-to-para ratio using the kinetic temperature obtained 
by the -R(2,2) /(i,i) • The ortho-to-para ratio depends on the 
temperature when NH3 molecules are formed. The ortho- 
to-para ratio increases in the case that NH3 molecules are 
produced in low temperature conditions (e.g., more than 
10 at 5 K) (Takano et al. 2002). i?( 3 , 3 )/(i,i) is the direct 
observed value that is related to the ortho-to-para ratio, 
and therefore, it should be studied in detail. 

The (3,3) emission is stronger than the (1,1) emission 
at most observed positions. The NH3 integrated inten- 
sity for i?(3 i 3)/(i.i) > 1 is 90% of the total NH3 intensity. 
The histogram of ii(3 3)/(i 1) is shown in Figure 5(b). The 
mean value and standard deviation of R(s,3)/n,i) are de- 
rived to be 1.28 and 0.29, respectively. R(3^)/(i^i) of the 
CMZ is higher than that of typical molecular clouds in the 
galactic disk. In galactic disk clouds, i?(3 ) 3)/(i,i) as high 
as that observed in the CMZ is found only in active and 
massive star- forming regions (see section 4.2). 

To investigate the spatial distributions of i?(2,2)/(i,i) 
and J2( 3i 3)/(i 1 i), we show the l-v diagrams of -R(2,2)/(i,i) 
and -R(3.3)/(i.i) superimposed on the (1,1) emission inte- 
grated over the entire observed latitude in Figure 6. The 
values of -R(2,2)/(i,i) and i?(3,3)/(i,i) vary with the loca- 
tion. The high-ratio region is common in both the ratios. 
Although the kinematical timescalc significantly depends 
on the galactocentric distance, neither ratio exhibits any 
systematic structure dependent on the distance from the 
Galactic center. Overall, the distributions of -R(2,2)/(i,i) 
and -R(3.3)/(i.i) are similar. However, some regions such 
as the Sgr C cloud near I ~ — 0.°5 exhibit a much lower 

#(3,3)/(l,l)- 

4. Discussion 

4-1. Physical Properties of Molecular Gas in CMZ 

Tk derived from the i?(2,2)/(i,i) IS thought to be the 
current temperature of the gas. The temperature derived 
from the ortho-to-para ratio corresponds to the tempera- 
ture when the NH3 molecules were formed (Umcmoto et 
al. 1999). These two temperatures can be obtained from 
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the NH3 observations. In the case that the entire interstel- 
lar matter is under thermal equilibrium and in a steady 
state, the two temperatures are the same. 

We derived the Tk of the gas under the cases of optically 
thin (t <C 1) and thick (r ~ 3) using the same method as 
that described in Paper I. The gas is warmer than the 
dust (Tdust = 15-22 K; Lis et al. 2001) at most observed 
positions, which is consistent with the conclusion shown 
in Paper I. The i?( 2 ,2)/(i,i) = 0.4-0.8 and i?( 2 ,2)/(i,i) > °-8 
correspond to Tk ~ 20-80 K and Tk > 80 K, respectively. 
The gases corresponding to Tk values of 20-80 K and 
> 80 K contain 82% and 18% of the total NH 3 integrated 
intensity in the entire observed area, respectively. The 
higher ratio (-R(2,2)/(i,i) > 0.7) components are located at 
I ~ -0.°5, 0.°0, 0.°9, and l.°3 (Figure 6(a)). These locations 
coincide with those of high CO J = 3-2/ J = 1-0 ratios 
(-^CO3-2/i-0 > 1-5; Oka et al. 2007). This supports the 
notion that the components with higher ii(2,2)/(i,i) values 
trace the high-temperature gas. -R(2,2)/(i,i) > 0.7 corre- 
sponds to Tk > 60 K in the optically thin case and > 40 K 
in the optically thick case (r ~ 3). This is consistent with 
i? CO 3-2/i-0 > 1-5 corresponding to Tk > 48 K, which is 
derived from a large velocity gradient (LVG) calculation 
(Oka et al. 2007). 

The high i?(3.3)/(i.i) value in the CMZ indicates that the 
ortho-NH3 is overabundant as compared to para-NH3. To 
obtain the ortho-to-para ratio, we estimated the column 
densities using the rotation diagram method under the op- 
tically thin assumption. We derived the column densities 
of each level from the integrated intensities. The proce- 
dure is based on Takano et al. (2002). Although the opti- 
cal depth of the NH3 lines in the CMZ is moderate (r < 
2-4; Hiittemeister ct al. 1993; Paper I), and the ortho- 
to-para ratio changes with the optical depth, the value 
obtained under the optically thin assumption can provide 
the lower limit of the ortho-to-para ratio (R / p = 2 at 
t<1 corresponds to R / p = 6 at r = 3). 

Because of the absence of emission lines from molecules 
at the (0,0) level, we cannot obtain the column density of 
the (0,0) level directly. However, we can derive the column 
density of the NH3 at the (0,0) level under the assump- 
tion that the ortho- and para-NH3 molecules are thcrmal- 
ized at the same rotational temperature derived from the 
-R(2.2)/(i,i)- This assumption is a better approximation 
than using the rotational temperature between (3,3) and 
(6,6), because the (1,1) and (2,2) levels are between the 
(0,0) and (3,3) levels. The energies of the (3,3) and (6,6) 
levels are 124.5 K and 412.4 K above the ground state, re- 
spectively. The large energy difference between the (3,3) 
and (6,6) levels suggests that the (3,3) and (6,6) emissions 
originate from different components. 

We calculated the column densities of ortho-NH3 and 
para-NH 3 by the formulas JV(ortho) = N(0, 0) + 2V(3, 3) 
and iV(para) = 2V(l,l)+iV(2,2), respectively. The column 
densities of the levels higher than (4,4) are less than 2% 
of those of the (0,0) and (1,1) levels when the rotational 
temperature is 36 K, which is estimated from the mean 
-R(2,2)/(i,i) = 0.70. We considered only the mctastablc 
states ( J = K), because transitions from non-metastable 



states (J 7^ K) to metasable states are generally quite 
rapid. 

Figure 7 shows histograms of the ortho-to-para ratio 
for the integrated intensity over the entire observed area 
with 10 a detection in the (1,1), (2,2), and (3,3) lines 
after smoothing to 10 km s _1 . The mean value and the 
standard deviation of R / p are derived to be 2.0 and 0.6, 
respectively. To investigate the spatial distributions of 
the ortho-to-para ratio, we have plotted the ortho-to-para 
ratio on the l-v plane, as shown in Figure 8. We found 
that the ortho-to-para ratio varies with the location, and 
it does not exhibit systematic structures dependent on the 
distance from the Galactic center nor the core/envelope of 
clouds. 

The distributions of i?(2,2)/(i,i) and i?(3,3)/(i,i) appear 
similar on the l-v plane (Figure 6). To see this correlation 
quantitatively, we made a correlation plot of -R(2,2)/(i,i) 
and -R(3,3)/(i,i), as shown in Figure 9. #(3 3)/(i,i) depends 
on both R /p and T rot , although -R(2,2)/(i,i) depends on 
only T rot . Therefore, we can draw curves of constant R / p 
values on the -R(2,2)/(i,i)"-^(3,3)/(i,i) domain (Figure 9). 
Figure 9 shows the curves with R / p = 1, 2, 4, and 6. 
Figure 9 shows also a plot of the observed ratios of the 
clouds in the galactic disk obtained from other observa- 
tions. Most galactic disk clouds such as M17, Cep A, and 
NGC 7538, show R / p ~ 1, although intense and massive 
starforming cores such as W49 and W51 are located along 
the curve with R a / p ~ 4. The CMZ is close to the latter 
and is different from the former. This suggests that the 
entire CMZ should be under physical conditions as similar 
to an intense starforming core. 

It has been reported that the ortho-to-para ratios of 
nearby external galaxies are as high as that of the CMZ. 
The ortho-to-para ratios of Maffei 2, Arp 220, and NGC 
253 arc 2.6, 1.7, and 6.2, respectively (Takano et al. 2000; 
2002; 2005). We calculated the orhto-to-para ratios of IC 
342 and M 82 using the data of Mauersberger et al. (2003) 
and Weifi et al. (2001), respectively. These have been plot- 
ted in Figure 9. From these data, IC 342 exhibits lower 
R / p as compared to that of other galaxies. However, 
i? / p of IC 342 is calculated to be 1.8 ± 0.5 from the data 
obtained using the NRO 45 m telescope (S. Takano, 2008 
private communication). Therefore, IC 342 may actually 
have a high R a / P value. 

M82 also shows low i? / P - Usero et al. (2007) suggests 
that the evolutionary stage of the starburst in M 82 is dif- 
ferent from those of NGC 253 and IC 342. This difference 
may relate to the low ortho-to-para ratio in M 82. 

The ortho-to-para ratio is > 2 for the CMZ and most 
of external galaxies. We conclude that the ortho-to-para 
ratio of the CMZ and the central region of the galaxies is 
higher than that of galactic disk. 

4.2. Origin of NH 3 

The ortho-to-para ratio is related to the temperature 
when NH3 molecules are formed. The ortho-to-para ra- 
tio approaches unity when NH3 molecules are produced 
under high temperature (> 40 K) conditions. This ratio 
increases when NH3 molecules are produced and equili- 
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brated at low temperature (e.g., more than 10 at 5 K). 
This is because most of molecules then reside in the lowest 
state, or the (0,0) ortho level. The relationship between 
the ortho-to-para ratio and the formation temperature is 
shown in Figure 3 of Takano et al. (2002). Our observa- 
tion shows the ortho-to-para ratio is 1.5-3.5 for 95% of 
the total integrated intensity of the entire CMZ (Figure 
7). This ratio corresponds to a formation temperature of 
11-20 K. This is much lower than the gas kinetic temper- 
ature derived from -R(2,2)/(i,i)- 

The low formation temperature of NH3 suggests that 
most of the NH3 in the CMZ was not produced by a gas- 
phase reaction. If NH3 were produced with a formation 
temperature of > 20 K, that is, the current gas tempera- 
ture, the ortho-to-para ratio should be < 1.5. Only 5% of 
the total integrated intensity shows R / p < 1.5. 

It should be noted that the formation temperature of 
NH3 is similar to the dust temperature in the CMZ, which 
is estimated to be ~ 15-22 K (Lis ct al. 2001). This 
suggests that the observed NH3 molecules were formed on 
the surface of cold dust grains. 

At the observed formation temperature, i.e., 11-20 K, 
NH3 molecules cannot be released from the dust surface 
into a gas phase. Therefore, the dust should be heated 
up to 100 K to release the NH 3 from the dust to the gas 
phase, because the ice at 100 K can release NH3 within 
a few years (Sandford & Allamandola 1993). The tran- 
sition timescale of ortho-NH3 and para-NH3 is ~ 10 6 yr 
in the interstellar gas (Cheung et al. 1969). Since the 
observed ortho-to-para ratio is larger than the statistical 
equilibrium value, NH3 would be released from the dust 
within the transition timescale. Therefore, the dust must 
be heated within < 10 6 yr. 

4-3. Transient Heating of Dust 

In the previous section, we suggested that the bulk of 
NH3 gas in the CMZ must be released from the dust that 
was previously heated up. We attribute this heating pro- 
cess to the followings possibilities: interstellar shocks, di- 
rect dust collisions, a starburst, and ambipolar diffution. 
Considering these mechanisms, we suggest that the inter- 
stellar shocks caused by a supernova and the direct dust 
collisions would be possible mechanisms for our transient 
dust heating scenario. 

We first consider the passage of an interstellar shock 
(Umemoto et al. 1999). For galactic disk objects, the en- 
hancement of the ortho-to-para ratio due to an interstellar 
shock is supported by the NH3 observations of the active 
and massive star- forming regions (Figure 9). In such a 
source, a strong bipolar outflow is produced, and it creates 
a strong shock, although the shocked region is much more 
compact than the CMZ. The interstellar shock can heat 
the interstellar gas directly. Although the gas can hardly 
heat the interstellar dust, accelerated electrons can heat 
the dust grain. 

In the CMZ, there is considerable evidence suggesting 
the presence of an extended shocked region. The ubiq- 
uity of the intense thermal SiO lines (Martin-Pintado et 
al. 1997; Handa et al. 2006) and many shell structures on 



the l-v map observed in the CS line (Tsuboi et al. 1999) 
are attributed to be caused by the strong shocks pass- 
ing through the entire CMZ. These shocks can transiently 
heat the dust. 

What is the origin of the shock? Multiple supernovae 
(SN) shocks are the first possibility. We estimated the 
timescale at which SN shocks blow through the CMZ to 
be ~ 10 4 yr using the multi-phase interstellar gas model 
of McKee & Ostrikcr (1977). In this estimation, we used 
an SN rate of lO^-lO" 9 pc" 3 yr" 1 (Muno et al. 2004; 
Schanne et al. 2006), and the H atom density at the SN 
remnant interior of — 10~ 2 cm -3 (McKee & Ostriker 
1977). In the Galactic center region, the H atom density 
of the ambient interstellar medium is larger than that in 
the galactic disk region, the H atom density at the SN 
remnant interior may also be larger. However, even if 
rih = 10 2 enr 3 , the timescale is ~ 10 5 yr. 

Galactic bar shock is the second possibility. The molec- 
ular gas in the CMZ should be strongly affected by a 
non-axisymmetric potential. The gas forms a bar-like 
structure, and is undergoing a strong non-circular motion 
(Sawada et al. 2004). In such a case, a strong shock is ex- 
pected at the leading edge of the bar (e.g., Sorensen et al. 
1976; Athanassoula 1992). This shock is observed in extra 
galaxies (e.g. Handa et al. 1990; Usero et al. 2006). We 
estimated the timescale at which the galactic bar shock 
passes through the entire CMZ to be 4 x 10 6 yr which is 
half of the rotation period at a radius of 250 pc and a 
velocity of 200 km s -1 . However, the absence of system- 
atic structure in the ortho-to-para ratio suggests that the 
galactic bar does not produce the fresh NH 3 . 

Cloud-cloud collisions are also a possible trigger of the 
shock wave. The energy of turbulent motions in the clouds 
within 300 pc of the center is 10 53 erg (Guesten et al. 
1985). This energy is larger than the thermal energy re- 
quired to heat the dust to 100 K, 10 51 erg, which is es- 
timated from the total dust mass in the central 400 pc 
of 5.3 x 10 7 M Q (Pierce-Price et al. 2000), and the spe- 
cific heat of graphite of 1.3 x 10 8 erg K _1 g~ 1 . Cloud- 
cloud collisions may heat the dust surface by the trans- 
fer of energy through the shock wave. However, the 
timescale of cloud-cloud collisions is one order of mag- 
nitude longer than the transition timescale of NH3, ~ 10 6 
yr, as shown below. The collision timescale is estimated 
by t co ] = ^/(n c o~cVc), where n c is the number density of 
clouds, cr c is the collisional cross section, and v c is the rel- 
ative velocity. The number density of clouds is estimated 
from the total molecular mass in the entire CMZ of (3- 
8) xlO 7 M Q (e.g. Tsuboi ct al. 1999) and a cloud mass 
of 2 x 10 3 Mq. We use a cloud radius of 1 pc and a H2 
number density of 10 4 cm' 3 in this estimation because 
clouds with a size of several pc are observed in the CS 
maps of Tsuboi et al. (1999), and the critical density of 
CS is 10 4 cm~ 3 . The collisional cross section is <r c ~ 10 37 
cm 2 at a cloud radius of 1 pc. The relative velocity of 
clouds is assumed to be v c ~ 10 km s _1 . The timescale 
of cloud-cloud collisions is estimated to be ~ 10 7 yr from 
these values. The comparison of the timescales suggests 
that cloud-cloud collisions cannot heat the dust frequently 
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enough. 

The direct collision of dust particles is the second mech- 
anism. The dust surface is heated to ~ 10 3 K when 
all of the collisional energy is transferred to the heat- 
ing energy of dust. In this mechanism, the dust particles 
must collide within the conversion timescale of ortho-NH 3 
and para-NH3. The collision timescale is estimated by 
icol = l/("dOdVd), where is the number density of dust 
particles, dd is the collisional cross section, and Vd is the 
turbulent velocity. The total dust mass within the central 
400 pc is 5.3 x 10 7 M Q (Pierce-Price et al. 2000), and the 
mass of a dust particle is estimated to be 6 x 10 -12 g for 
a dust particle radius of 1 /zm and mass density of char- 
coal powder of 1.5 g cm -3 . Therefore, the number density 
is estimated to be rid ~ 10~ 10 cm~ 3 . This is an average 
value over the entire CMZ. The collisional cross section is 
o"d ~ 10~ 8 cm 2 at a radius of 1 /zm. The turbulent velocity 
is assumed to be Vd ~ 10 km s _1 . Using these values, we 
estimated t co \ ~ 10 4 yr. This timescale is upper limit. If 
the dust is concentrated in the core of each cloud with the 
same turbulent motion, rid is larger and it gives shorter 
t co i. Before the grains collide directly, the gas clouds that 
contain dust particles collide with each other. In this case, 
the dust particles can slow down before the direct collision. 
However they are also affected by gas drag, which heats up 
the dust. Therefore, the conclusion remains valid, unless 
the dust particles are tightly bounded in the molecular 
clouds by the gas drag. In this case, the discussion should 
be the same as that of cloud-cloud collisions. The com- 
parison of timescales suggests that the direct collision of 
dust particles is a possible heating mechanism. 

Another heating mechanism is a transient starburst. 
When a dust particle is heated up to 100 K by strong 
UV radiation, NH3 molecules are released in the interstel- 
lar space. There is substantial evidence that indicates a 
burst of star formation in the Galactic center. However, 
this occurred 10 7 -10 9 years ago (Yasuda et al. 2008 and 
reference therein). This timescale is much longer than 
the transition timescale of NH3. Therefore, a transient 
starburst could not have heated the dust. 

The last mechanism is an ambipolar diffusion. 
Hiittcmcistcr et al. (1993) suggested that the gas can be 
heated up to ~ 200 K by collisions between neutrals and 
ions under a magnetic field strength of 500 pG. However, 
this mechanism does not work well in our transient heat- 
ing scenario. The energy of collision between the dust 
particle and an ion is fifteen order of magnitude smaller 
than that of direct dust collision. The collisional energy 
is too small to heat the dust to a sufficient temperature 
to release the NH3, 100 K. Therefore, the dust could not 
been heated by ambipolar diffusion, 

Our model requires the heating to have occurred within 
< 10 6 yr. Therefore, SN shocks, galactic shocks, and direct 
dust collisions are possible. However, our estimation of the 
cooling timescale of molecular gas, which is described in 
the next subsection, indicates that the timescale should 
be ~ 10 5 yr. Therefore, a galactic shock that occurred 
~ 10 6 yr ago is not preferable. We therefore conclude 
that SN shocks and direct dust collisions would be possible 



mechanisms for our transient dust heating scenario, and 
both mechanisms can produce fresh NH3 from the dust 
grains. 

We expect that there should be the difference of the 
spatial distributions in the SN shocks and the direct dust 
collisions. In the case of the SN shocks, the dust would be 
heated in the shocked layer. In the case of the direct dust 
collisions, individual heating spot is the fairly smaller scale 
of dust particles. These minor heated dusts are mixed 
with the major cool dusts in the observed beam, therefore 
they arc uniformly observed as cold dusts. Submillimctcr 
wave observations with very high resolution could find the 
differences of the distributions between two mechanisms. 

4-4- Cooling Timescales of Gas and Dust 

NH3 is released from the dust grains to the gas-phase 
at a sublimation temperature of ~ 100 K. Therefore, the 
dust should be heated up to to ~ 100 K in the case that 
NH3 molecules are released from the dust. The observed 
gas temperature is close to this temperature. However, 
the observed dust temperature is much lower. 

To resolve this discrepancy, we hypothesize that the in- 
terstellar matter in the Galactic center region is not in 
thermal equilibrium. The dust is believed to heat the gas 
through the photoelectric process, however, the gas can- 
not heat the dust. Therefore, our hypothesis would be 
confirmed if the cooling timescale of gas is much longer 
than that of dust. 

We estimated the timescale in which the gas is cooled 
from 100 to 10 K. The main cooling process of the gas 
in this temperature range is CO emission. The cooling 
timescale of gas is estimated to be ~ 10 5 yr from the ther- 
mal energy of Q — 5 /2k AT ~ 10 -14 erg and the CO cooling 
rate of ~ 10~ 27 erg s" 1 (Hegmann et al. 2007). 

We estimated the timescale in which the dust is cooled 
from 100 to 10 K. The major cooling process of the dust 
in this temperature range is black body radiation. The 
thermal energy of the dust with temperature T(t) at time 
t is given by 

Air 

Q(t) = —a 3 pcT(t) (1) 

where a, p, and c denote the radius, density, and specific 
heat of dust. The cooling rate of the dust by the black 
body radiation is given by 

dQ(t) 



dt 



= -4ira 2 a[T(t)Y 



(2) 



where a is the Stefan-Boltzmann constant. Using these 
equations, we found the cooling timescale of dust from 
100 to 10 K to be t ~ 10" 3 yr for a = l^m, p = 1.5 g cm -3 
as the mass density of charcoal powder, and c = 1.3 x 10 8 
erg K~ x g" 1 as the specific heat of the graphite. 

The cooling timescale of the gas is 7 orders of magnitude 
longer than that of the dust. Therefore, transient heating 
of the interstellar matter occurred in the past, and the 
observed temperature difference between the gas and the 
dust can be observed soon after the heating. This may 
explain why temperature of molecular gas is much higher 
than that of the dust. 
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4-5. Individual Clouds 

4-5.1. Overview of Clouds in the CMZ 

We estimated the physical parameters of the identi- 
fied 13 clouds under the optically thin assumption. The 
calculated rotational temperatures, column densities, and 
ortho-to-para ratios are listed in Table 3. 

With regard to the clouds in the CMZ, the virial masses 
observed in the CO line are one order of magnitude larger 
than the mass derived from the CO luminosity with a 
standard conversion factor (Oka et al. 1998). This sug- 
gests that the clouds in the CMZ may not be gravitation- 
ally bound and are instead in pressure equilibrium with 
the hot gas or magnetic field in this region. However, the 
cloud observed in the NH 3 line may not be affected by 
the external pressure because the NH3 line would trace a 
different density range. 

Therefore, we investigate the virial mass and luminos- 
ity mass in the NH3 line. Using the abundance ratio of 
X(NH 3 ) = 10~ 9 (Huttcmeister ct al. 1993) and spherically 
symmetric geometry, the H 2 number densities are derived 
from the luminosity masses to be n(H 2 ) ~ 10 3 -10 4 cm~ 3 , 
as given in Table 3. These values are consistent with 
the NH3 critical density, suggesting that the estimated 
luminosity masses are valid. We estimated the luminosity 
mass, M\ um , from the integration of the column densities, 
and the virial mass, M v j r , from the apparent size with 
the beam deconvolution and the observed line width. We 
summarize these values of the 13 clouds in Table 3 and 
show the M\ um -M vir plot in Figure 10. The luminosity 
mass and virial mass of nine clouds are consistent with an 
order of magnitude. This implies that most of the molec- 
ular clouds traced by the NH3 line are roughly in virial 
equilibrium. The external pressure affects only the less 
dense gas traced by the CO line. 

Only four clouds (ID 1, 3, 6, and 8) appear to not be in 
the virial equilibrium, although there are reasons in each 
case, as given below. In the l-v diagram, these clouds do 
not lie on the main ridge that indicates the main rotating 
component of the Milky Way Galaxy. The rotational tem- 
perature and the ortho-to-para ratio of these clouds are 
different from those of the other clouds. Wc summarize 
the features of these four unvirialized clouds as follows. 
4.5.2. Cloud 1 

Cloud 1 is located at the western edge of the observed 
area and at «lsr — km s _1 . The velocity and angular 
extent suggest that cloud 1 may be located in the Galactic 
disk. The annihilation source 1E1740. 7-2942 is located 
close to the cloud in the sky. Bally & Leventhal (1991) 
assumed that the source is not associated to the km s _1 
component because it should not be in the Galactic center 
region. However, the rotational temperature of 45 ± 4 K 
is higher than that of the disk clouds and close to that of 
the clouds in the Galactic center. Therefore, we conclude 
that cloud 1 may be located at the Galactic center region. 
In a higher-resolution map in the CS line (Tsuboi et al. 
1999), two clouds are observed at the position of cloud 1. 
In the case that these clouds are gravitationally unbound, 
our estimation of the virial mass is overestimated and the 



large discrepancy between the two masses is not real. 

4.5.3. Cloud 3 

Cloud 3 is located at Sgr C in the sky. The peak ve- 
locity in the NH 3 line of i>lsr — —53 km s _1 is close 
to the central velocity in the H91a and H109a recom- 
bination lines of «lsr — —59 km s _1 (Pauls & Mezger 
1975; Anantharamaiah & Yusef-Zadeh 1989). This sug- 
gests that the Sgr C Hll region is physically associated. 
It interacts with the molecular cloud and ionizes a part of 
it. The rotational temperature of 45 ± 3 K is higher than 
that of other clouds in the CMZ. This suggests that cloud 
3 is heated by the H II region and it is expanding rapidly 
due to the active star formation. 

4.5.4. Cloud 6 

Cloud 6 is located at high Galactic latitude (b ~ 0.15) 
at ^lsr ~ 106 km s _1 . It is redshifted approximately 60 
km s _1 from the Sgr A 50 km s _1 cloud, which is located 
at the same longitude. This cloud appears as a shell-like 
structure on the l-v map in the CS line (Tsuboi et al. 
1999). Soft X-ray (1-3 keV) emission is detected at the 
position of this cloud (Wang et al. 2002). These suggest 
that this cloud is associated with a newly identified super- 
nova remnant. The large virial mass should be due to ex- 
pansion of this SN remnant. The kinetic energy is 10 51 erg 
when cloud 6 is expanding with a mass of M\ um = 2.1 x 10 5 
M Q and an expansion velocity of Aw/2 = 17 km s _1 . This 
energy is comparable to the kinetic energy of the stellar 
wind from a 60-100 M Q star (Abbott 1982) and less than 
that of a hypernova explosion (Iwamoto et al. 1998). 

4.5.5. Cloud 8 

Cloud 8 is the same structure reported in Paper I. The 
rotational temperature and the number density of this 
cloud are higher than those of other clouds. The location 
of this cloud is close to that of the 0.°9 anomaly reported 
from the observation in CO line (Oka et al. 2007). The 
features of high temperature and density are similar to 
the 0.°9 anomaly. Although the highly blue-shifted com- 
ponent (wlsr = —120 km s _1 ) is detected in CO line, we 
could not detect this component. 

This cloud is unresolved by the Kagoshima 6 m tele- 
scope beam. However, it could be resolved by a high res- 
olution observation using the Kashima 34 m telescope. 
Figures 11 and 12 show the l-v and b-v maps obtained us- 
ing the Kashima 34 m telescope. In this higher-resolution 
observation, the cloud is resolved into three components. 
The peak positions, apparent sizes, and line widths of 
these components arc shown in Table 4. 

We call these three components 8a, 8b, and 8c in the 
order of the galactic longitude. Cloud 8 has a large virial 
mass because it is an unbound system. The rotational 
temperatures of all components are similar, and they are 
higher than the mean temperature of the CMZ. The lo- 
cation and velocity of components 8c is close to the high 
velocity compact cloud, CO 0.88-0.08 (Oka et al. 2007). 
These would be the same structure. However, the com- 
ponent 8a and 8b are not detected in CO line (Oka et 
al. 2007). This would be because the CO line is affected 
by foreground absorption. We summarize the physical pa- 
rameters of three components that have been obtained un- 



8 



T. Nagayama et al. 



der the optically thin assumption in Table 5. Their virial 
masses are the same as their luminosity masses within a 
factor of ~ 2. 

5. Conclusions 

We have presented a map of the major part of the CMZ 
by simultaneous observations in the NH3 (J, A') = (1,1)) 
(2,2), and (3,3) lines using the Kagoshima 6 m telescope. 
Our observations are summarized as follows: 

1. Wc obtained the ortho-to-para ratio of NH3 in 
the CMZ to be 1.5-3.5 using the rotation diagram 
method. This ratio is higher than that of galac- 
tic disk clouds, and as same as those of the active 
and massive star forming regions and the centers of 
nearby external galaxies. 

2. An ortho-to-para ratio of 1.5-3.5 corresponds to a 
formation temperature of 11-20 K. The formation 
temperature is very close to the temperature of dust. 
NH3 in the CMZ has been produced on the cold dust 
grains and released into the gas phase by SN shocks 
or direct dust collisions. 

3. The discrepancy between warm molecular gas and 
cold dust can be explained by transient heating of 
the interstellar media in the CMZ on a timescale 10 5 
years. NH3 molecules were released from the dust 
on this timescale. 

4. We present the physical conditions of the 13 clouds 
identified in the NH3 map. For a majority of the 
clouds in the CMZ, the luminosity mass and the 
virial mass are consistent in an order of magnitude. 
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Table 1. Molecular clouds toward the Galactic center identified in the NH3 line. 



ID 


I b 


^LSR 


Al 


Ab 


Av 


Comments 




(°) (°) 


(km s _1 ) 


(°) 


(°) 


(km s _1 ) 




1 


-0.86 -0.11 


-10 




0.21 


53 


at the western edge of the observed area 


2 


—0.86 0.00 


141 


0.20 


0.17 


31 


a part ol .LMR 


3 


-0.51 -0.15 


-53 


0.17 


0.16 


39 


associated with Sgr C 


4 


-0.11 -0.07 


21 


0.21 


0.16 


49 


Sgr A 20 km s _1 cloud 


5 


0.07 -0.08 


48 


0.23 


0.15 


49 


Sgr A 40 km s" 1 cloud 


6 


0.07 0.15 


106 


0.16 


0.19 


33 




7 


0.74 -0.13 


39 


0.34 


0.27 


65 


Sgr B cloud complex 


8 


0.78 -0.12 


22 




0.16 


56 


0.°9 wing feature* 


9 


0.97 -0.10 


80 


0.31 


0.19 


37 




10 


1.24 -0.07 


83 


0.43 


0.35 


49 


the l.°3 cloud complex 


11 


1.53 -0.23 


-25 


0.23 


0.17 


37 




12 


1.60 -0.07 


50 




0.24 


29 




13 


1.70 -0.37 


-35 


0.18 




27 





* defined in Paper I 



Table 2. The integrated intensity of 13 clouds. 



ID 


J TuBdv 


fraction* 




(K km s" 1 ) 


(%) 


1 


90 


1.8 


2 


12 


0.2 


3 


59 


1.2 


4 


305 


6.1 


5 


434 


8.7 


6 


58 


6.1 


7 


1224 


24.6 


8 


107 


2.2 


9 


246 


5.0 


10 


562 


11.3 


11 


135 


2.7 


12 


151 


3.1 


13 


64 


1.3 


total 


3447 


69.4 



* The traction ot the total integrated 
intensity of the entire observed area 
( / T MB dv = 4972 K km s" 1 ). 



NH 3 emission in CMZ 



11 



Table 3. Physical properties of the 13 clouds. 



ID 


T ro t 


iV(ortho) 


TV(para) 


AT(total) 




Ro/p 


-'form 


d 


n(H 2 ) 










(K) 




10 14 (cm" 2 ) 






(K) 


(PC) 


10 4 (cm" 


- 3 ) 


10 6 (M ) 


10 6 (M ) 


1 


45 ±3 


1.8±0.2 


1.69 ±0.09 


3.5 ±0.3 


1 


,1±0.2 


> 24 


21 


0. 


.55 


1.2 


7.2 


2 


31 ±2 


2.9±0.6 


1.08 ±0.07 


3.9±0.6 


2 


.6 ±0.7 


13 ±3 


13 





.98 


0.56 


1.6 
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Table 4. Three components of the 0.°9 wing feature observed using Kashima 34 m telescope. 



ID 


I 


b 


^LSR 


Al 


Ab 


Av 




(°) 


(°) 


(km s" 1 ) 


(°) 


(°) 


(km s" 1 ) 


8a 


0.767 


-0.030 


25 


0.065 


0.065 


28 


8b 


0.775 


-0.098 


22 


0.077 


0.090 


29 


8c 


0.854 


-0.059 


10 


0.061 


0.067 


35 



Table 5. Physical properties of the three components of the 0.°9 wing feature. 



ID 


?rot 

(K) 


TV(ortho) 


iV(para) 
10 14 (cm" 2 ) 


TV(total) 


R /p 


-'form 

(K) 


d 
(pc) 


n(H a ) 
10 4 (cm- 3 ) 


M\um 

10 6 (M Q ) 


10 6 (M Q ) 


8a 


40 ±1 


14.7±0.6 


9.98±0.18 


24.7 ±0.8 


1.9±0.1 


16 ±1 


9 


9.1 
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0.86 


8b 


40 ±1 


15.4±0.6 


10.7±0.18 


26.1 ±0.8 


1.5 ±0.1 


20 ±2 


12 


7.3 


3.0 


1.2 


8c 


37 ±1 


19.6 ± 1.1 


9.11 ±0.21 


28.7±1.3 


2.2±0.2 


14 ±1 


9 


11 


1.8 


1.3 
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Fig. 1. The NH3 (1,1) (red), (2,2) (green), and (3,3) (blue) spectra toward four positions of the observed area, Sgr C 
(I, b) = (-0.°500, -0.°125), Sgr A (0.°125, -0.°125), Sgr B (0.°750, -0.°125), and l.°3 region (l.°250, -0.°125). Spectra are 
shown with the veolcity resolution of 5 km s -1 and on a Tmb scale. 
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Fig. 2. Integrated intensity maps of NH3 (1,1) (left), (2,2) (middle), and (3,3) (right). The lowest contour and the contour interval 
are 2,2 and 4.4 K km s — 1 in I Tyi^dv, respectively. The velocity range of the integration is —200 < ulsr 5: 200 km s . 



200 




/(degree) /(degree) /(degree) 



Fig. 3. Longitude-velocity (l-v) diagrams of NH3 (1,1) (left), (2,2) (middle), and (3,3) (right) emission integrated over the entire 
observed latitude. The lowest contour and the contour interval are 0.25 and 0.31 K, respectively. 
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Fig. 4. Longitude-velocity (l-v) diagrams of the NH3 (1,1) (left), (2,2) (middle), and (3,3) (right) lines. Both of the lowest contour 
and the contour interval are 0.1 K. 
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Fig. 4. Continued 
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Fig. 5. Histograms of the intensity ratios of (2,2) to (1,1) (left), and (3,3) to (1,1) (right). We count the l-b-v pixels at which the 
lines were detected above the 3 <r level after smoothing at 10 km s _1 . 
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Fig. 6. The distributions of H(2,2)/(l,l) ( a ) an d ^(3,3)/(i,l) (b) superimposed on the NH3 (1,1) emission integrated over the entire 
observed latitude. T ro t in (a) is derived for the optically thin case. T ro t in (b) is derived for the optically thin case and an 
ortho-to-para ratio of 2.0. 




Fig. 7. The histogram of the ortho-to-para ratio. We count the l-b-v pixels at which the lines were detected above the 10 a level 
after smoothing at 10 km s . 
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Fig. 8. The distributions of R / p superimposed on the NH3 (1,1) emission integrated over the entire observed latitude. 
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Fig. 9. The R(a 2)/(i i)~R(s,3)/(l,l) correlation plot of the CMZ (black), galactic disk clouds (green), and the central region of 
some external galaxies (red). For the CMZ, the spectral data detected over the 10 <r level after smoothing at 10 km s — 1 are plotted. 
The data of galactic disk clouds and external galaxies are referred to from the following; Cep A (Ho et al. 1982); M17 N, M17 NH3, 
and M17 S (Guesten & Fiebig 1988); NGC 2264 and NGC 7538 (Nagayama et al. 2008 in prep.); Orion KL (Ho et al. 1979); S140 
(Mauersberger et al. 1985); W3 West and W3 OH (Tieftrunk et al. 1998); W49 (Nagayama et al. 2008 in prep.); W51 (Matsakis et 
al. 1980); W75 Field 1, W75 NH 3 -5, and W75 S (Wilson k, Mauersberger 1990); Arp 220 (Takano et al. 2005); IC 342 (Mauersberger 
ct al. 2003); Maffci 2 (Takano et al. 2000); M82 (Weifi et al. 2001); NGC 253 (Takano et al. 2002); Arp 220 is plotted from the ratio 
of optical depths in the absorption lines. 
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Fig. 10. The luminosity mass-virial mass plot for the 13 clouds. 
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Fig. 11. High-resolution l-v diagrams of cloud 8 observed using the Kashima 34 m telescope in the NH3 (1,1) (left), (2,2) (middle), 
and (3,3) (right) lines. Both the lowest contour and the contour interval are 0.2 K. 
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Fig. 12. High-resolution b-v diagrams of cloud 8 observed using the Kashima 34 m telescope in the NH3 (1,1) (left), (2,2) (middle), 
and (3,3) (right) lines. Both the lowest contour and the contour interval are 0.2 K. 



